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Abstract 

We propose a theoretical framework to describe the ladder systems. The N-chain 
Hubbard model has been studied within the Composite Operator Method. In this 
scheme of calculations the single-particle Green's function for any number of coupled 
chains is obtained by solving self-consistently a system of integral equations. 
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The study of the so-called ladder materials is an actual hot field in Con- 
densed Matter Physics. In these systems quantum effects become very im- 
portant and lead to a dramatic dependence on the number of coupled chains 
that constitute the ladder [1]. Physical realizations of such systems are, for 
instance, vanadyl pyrophosphate ((VO)2P207) and the family of layer com- 
pounds Sr„_iCu n+ i02n- Even-leg ladders have a gap in the spin excitation 
spectrum whereas odd-leg ladders present no spin gap [1]. Such difference is 
confirmed by measurements of the magnetic susceptibility and neutron and 
muon spin scattering [2]. Ladder systems are studied by different models like 
the t-t'-J-J' [3], Heisenberg [4] or Hubbard [5] models. In this paper we propose 
a theoretical framework for these systems. We consider N coupled Hubbard 
chains and solve the model in the Composite Operator Method (COM). 

We consider a two-dimensional lattice described by the lattice vectors 

Ri = xi x + yi y i x = na [— oo < n < oo] i y — mb [1 < m < N] (1) 



For open boundary conditions along the ^-direction, the Hubbard model is 
described by the following Hamiltonian 

# = E (tiJ ~ ct «c(j) + X>T(0»i(0 (2) 

ij i 



where c(i) denotes the electron field at the site Rj in a spinor notation; n a (i) = 
cJ.(i)c CT (i) is the density operator for electrons with spin a; U is the on-site 
Hubbard interaction and /i is the chemical potential. The hopping matrix Uj 
is given by i - = -2t x af xj J iyjy - 2t y a y iyjy S ixjx where t x and t y are the hopping 
amplitudes along and between the chains, respectively. When iV = 1 we put 
t y = 0, in order to recover the 1-chain model. cx% x j x and aOl j are the projection 
operators along the x- and y-direction, respectively: 
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For the prescribed boundary conditions the electron field has the following 
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representation 
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The Fourier transform of takes the expression 
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where t(k x , k y ) = —2t x cos(k x a) —2t y cos(k y b); a and b are the lattice constants 
along the x and y direction, respectively. In the framework of the COM [6] 
we consider the Hubbard doublet = (c'(i)(l — n(i)) , c* (i)n(i)) as the 

basic field, where n(i) = c\i)c(i) is the density operator. By means of the 
Hamiltonian (2) the Heisenberg equation for the composite field is given by 
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with c\i) = 'E j t ij c(j),n li (i) = c 1 (i) ■ a ^(i) , = ^{1)^(1) + c(i)[<* (i) 
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Let us consider the thermal retarded Green's function S(i,j) = (R[ip(i)ip^ (j)]) . 
In the two-pole approximation [6] we have the equation 



S(l,j)=i5(t i -t j )I(i,j) 



(9) 



where = ({^(i), ^(j)}) and m(i, 3) = V f 0')}) = e(M)A^)- 

By introducing the Fourier transform 
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Eq. (9) takes the following form in momentum space 

US(k x ,ky,Py,U) + Jf-p£2^e(k x ,ky,q y )S(k x ,q y ,Py,U) = l(ky,Py) (ll) 



ly 



The solution of this integral equation will determine the Green's function, once 
the matrices I(k y ,p y ) and e(k x , k y , q y ) are known. Lengthy but straightforward 
calculations show that the normalization matrix I{k y ,p y ) and energy matrix 
e(k x , k y , q y ) are determined by 



H k y,Py) = Y, sin ( k y i y sin (Py i y) I ( i y) 



(12) 
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e(k x , k y , q y ) = — — ]T m(k x , k y , q y ) ]T sm(q y i y ) sva{p y i y )r l {i y ) (13) 

iV + 1 Qy iy 



where 
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The elements of the m-matrix are given by 

m n (k X ,ky,Py) = -IJ,I n (ky,Py) + A(ky , Py) +t(k x ,ky) [hl(ky,Py) 

-I 2 2(ky,Py)] + P(k X , ky,Py) 
m i2(k X , ky, Py) = A(/Cj,, Py) + t(k X , ky) P(k X , ky , Py) 

m 22 (k X , ky,Py) = "(/I - U)l 22 (ky,Py) + A(ky , Py) + P (k x , ky , Py) 

with the following definitions 



(15) 



A ( k y,Py) =52fan(kyiy)sm(p y i y ) - rf(i)rf(i)) 

ly 

(21) P(k x ,ky,p y ) = -t x cos(k x a) \P x+ (k y ,p y ) + P x ~(k y ,p y ) 

~ty [py + {ky,Py)+Py-{k y ,P y ) 
P X± (ky, Py ) =Y,^(k y iy)s : m{j>yiy)P X± {iy) 
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P y± (ky, P y) =Y,Sm(kyly) SU1 i^yly ^ P y 6) (ly ) 

px± (h) = tKOK^Ta,^)) - (c](i)ci(i)c\(i x Ta,i y )c\(i x Ta,i y )) 
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Fig. 1. Chemical potential n as a, function of the ratio r = t y /t x for U = 4, T = 0.01 
and n = 0.9 in the two-leg ladder (N = 2). 

P v± (i y ) = \{n»{i)n^{i x ,i y ^b,iy)) 

-{^(i)ci(i)c\(i x , i y =F a, i y )c\(i x , i y T b, i y )) (18) 

As we see a series of unknown parameters appear. This is a consequence of 
the fact that the properties of the composite fields are not known a priori. 
They are determined by the dynamics and the boundary conditions and must 
be self-consistently calculated. Parameters such as /i, (n(i y )) and A(i y ) can be 
expressed in terms of matrix elements of the Green's function. The parameters 
P x (i y ) and P y± (i y ) are static spin, charge and pair correlation functions and 
can be determined by the content of the Pauli principle. We thus can write 
a series of coupled equations which lead to a self-consistent calculation of the 
Green's function. 

As preliminary result, the chemical potential as a function of the ratio r = 
t y /t x for the two-leg ladder system (N = 2) is presented in Fig. 1. The sharp 
feature around r 0.55 could indicate the onset of some spin ordering re- 
lated to the increasing of frustration in the spin coupling channel due to the 
appearance of an additional exchange interaction along the rungs. A deeper 
comprehension of the on-going dynamics requires the analysis of the correla- 
tion functions of the system. 
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